Distance-Redshift in Inhomogeneous Qq = 1 
Priedmann-Lemaitre-Robertson- Walker Cosmology 



R. Kantowski 

University of Oklahoma, Department of Physics and Astronomy, 
Norman, OK 73019, USA 

<^> kantowski@mail.nhn.ou.edu 

o 

^ : R. C. Thomas 

| ^ . University of Oklahoma, Department of Physics and Astronomy, 

Norman, OK 73019, USA 

thomas@mail . nhn . ou . edu 

(N 
> 

iP ! ABSTRACT 



o 

Distance-redshift relations are given in terms of associated Legendre functions 
for partially filled beam observations in spatially flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmologies. These models are dynamically pressure-free, flat FLRW 
on large scales but, due to mass inhomogeneities, differ in their optical properties. The 
partially filled beam area-redshift equation is a Lame' equation for arbitrary FLRW 
and is shown to simplify to the associated Legendre equation for the spatially flat, i.e., 
f2o = 1 case. We fit these new analytic Hubble curves to recent supernovae (SNe) data 
^ . in an attempt to determine both the mass parameter Q m and the beam filling parameter 

v. We find that current data are inadequate to limit v. However, we are able to esti- 
mate what limits are possible when the number of observed SNe is increased by factor 
of 10 or 100, sample sizes achievable in the near future with the proposed SuperNova 
Acceleration Probe satellite. 



Subject headings: cosmology: theory - large-scale structure of universe 



1. INTRODUCTION 

Distance-redshift or equivalently the Hubble curve is critical in determining current values of 
the cosmological parameters Hq, f2 m , and Oa- Conversely, current values of these three parameters 
determine the large scale dynamics of the Universe into the distant past. A complication occurs 
when attempting to determine these parameters from high z comparisons to the standard Hubble 
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curve. The standard Hubble curve is a theoretical quantity computed assuming all gravitating mat- 
ter is homogeneously distributed; whereas, observational data is taken in the real inhomogeneous 
Universe. In an inhomogeneous universe an observing light beam is lensed by inhomogeneities lo- 
cated external to, but near the light beam, and defocused (relative to the standard Hubble curve) 
by the less than average matter density within the beam. The simplest way to take into account 
these effects is to correct all beams for the missing homogeneous matter but correct for lensing only 
when necessary. This procedure requires the introduction of one additional parameter, e.g., a filling 
parameter v, < v < 2 defined by the fraction of inhomogeneous matter pi/ po = v{y + l)/6 < 1 
excluded from observing beams {y = is the standard 100% filled beam FLRW case and v = 2 
is the empty beam case). When observing high z objects [z ~ 1) the reader can think of the 
parameter v as representing matter that exists in galaxies but not in the intergalactic medium. To 
find the theoretical Hubble curve for observations in such a universe one must solve the geometrical 
optics equation [see Kantowski (1998)] given as equation (1) in the next section. This equation is 
actually equivalent to the Lame' equation for general FLRW but as pointed out by Kantowski et al. 
(2000) reduces to the associated Legendre equation (7) for the special case considered here, Qq = 1. 
In § 2 we solve this equation using appropriate boundary conditions and give the Hubble curve in 
terms of associated Legendre functions (eq.[10]) as well as in terms of hyper geometric functions 
(eq.[16]). In § 3 we fit this new Hubble curve to data for 60 supernovae (SNe) from the Supernova 
Cosmological Project (SCP) and from the Cala'n/Tololo Supernova Survey (CTSS) in an attempt 
to determine the mass parameter Q, m and the filling parameter v. In § 4 we give some concluding 
remarks. 



2. The Luminosity Distance-redshift Relation 

For models being discussed here (and for most cosmological models), angular or apparent size 
distance is related to luminosity distance by D < (z) = De(z)/(1 + z) 2 . We choose to give luminosity 
distances in this paper. The D((z) which accounts for a partially depleted mass density in the 
observing beam but neglects lensing by external masses is found by integrating the second order 
differential equation for the cross sectional area A(z) of an observing beam from source {z = z s ) to 
observer (z = 0), see Kantowski (1998) for some history of this equation, 1 : 

(1 + Z) Vl + tomZ + Ha[(1 + Z)' 2 ~ 1] X 

^(i + *)Vi + n mZ + n A [(i + z)-2 - 1] ^y/Mz) 

+ (3 + i/) 4 (2 ~"W + zf = 0. (1) 



1 This equation follows from applying Sach's optics equations [Sachs (1961)] to an inhomogeneous FLRW universe 
and neglecting shear (external lensing) Kantowski (1969). The first version of equation (1) was given by Zel'dovich 
(1964) and later Dyer & Roeder (1974) included the cosmological term. 
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The required boundary conditions are 
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where <5f2 is the solid angle of the beam at the source and the FLRW value of the Hubble parameter 
at z s is related to the current value Hq at z = by 



H s = H (l + z s ) ^/l + n m z s + n A [{l + z s )- 2 -l]. (3) 
The luminosity distance is then simply related to the area A\ Q of the beam at the observer by 



D; 



= 4ki 



(1+ZsY 



(4) 



Equation (1) can be put into the form of a Heun equation and its solution has been given in terms 
of Heun functions in Kantowski (1998). Even though the Heun equation is only slightly more 
complicated than the hyper geometric equation, e.g., it has 4 regular singular points rather than 
3, Heun functions are not yet available in standard libraries. Consequently, such expressions are 
not particularly useful for comparison with data, at this time. Because the exponents of three of 
the singular points of the area equation (in standard Heun form) are and 1/2 [see eq. (13) in 
Kantowski (1998)], equation (1) is actually equivalent to the doubly periodic Lame' equation. We 
now show that it reduces to the associated Legendre equation for the spatially flat universe, f2o = 1. 
The required change of dependent and independent variables are respectively 



P(A,z) EE (1+ Z f/^—, 



r](z) 



1 1 + n m z(3 + 3z + z 2 ) 
tt m (l + z)3 



(5) 
(6) 



The resulting associated Legendre equation is 



(1 - v 2 ) 



with initial conditions: 
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P\s = 
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The resulting luminosity distance is then given by 



D t (z a ) = (l + z s )P( V (0)). 



(9) 
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Expressed as associated Legendre functions equation (9) becomes 

c 2 r(^)r(^)(i + z) 3 / 4 



De(Q m , ft a = 

p (l+2z/)/6 
r -1/6 

(1+2j/)/6 
-1/6 
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(10) 



In this expression the associated Legendre functions take on their their analytically continued values 
(i.e., the arguments are on the real axis and > 1). When the filling parameter has values v = 0, 1, 
or 2, equation (10) reduces respectively to equations (22), (39) and (54), of Kantowski et al. (2000). 
Because the associated Legendre equation is a special type of the hypergeometric equation it is 
always possible to write associated Legendre functions in terms of hypergeometric functions. And 
because hypergeometric functions are the more universally available, these results are the more 
useful for most parameter values. That the hypergeometric result existed has independently been 
seen by Kantowski et al. (2000) and Damianski et al. (2000). That the area equation reduces to 
the associated Legendre equation for another special case (A = 0) has been known for some time, 
see Kantowski et al. (1995) and Seitz h Schneider (1994). The appropriate change of variables is 



h(A,z) = (1 + z) 
C(z) = 




= (l + z)-^P, 
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The resulting hypergeometric equation is 



d 2 h 



with initial conditions 
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The resulting luminosity distance is then given by 

D e (z s ) = (1 + z s )MC(0)). 
Expressed in terms of hypergeometric functions equation (15) becomes 

c (1 + z) 2 



(14) 



(15) 



D e (n m , n A = i - n m , w, z) 
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3. Prospects for constraining v from high redshift SNe la 

Recent cosmic microwave background observations strongly imply a spatially flat universe (de 
Bernardis et al. (2000), Roos & Harun-or-Rashid (2000)), which naturally motivates application of 
the new distance formulae presented in §2 to a set of standard candles to estimate v. We use the 60 
SNe la from the combined Calan/Tololo + Supernova Cosmology Project (CT+SCP) as presented 
in Riess et al. (1998), Perlmutter et al. (1999). Combining these data with those from the High-z 
SN search (Schmidt et al. (1998)) would bring the total to about 100 SNe. However, we shall see 
that this somewhat complicated task would not increase the numbers of SNe enough to noticeably 
improve the estimate. 

Rather than subject the data to an in-depth Bayesian re-analysis with the additional beam 
filling parameter v included, we merely use a x 2 goodness-of-fit estimation. We assume an intrinsic 
SN absolute magnitude of Mb = —19.33 and Hq = 65 km s _1 Mpc -1 and follow the same procedure 
as Wang (2000a) to recover results consistent with Perlmutter et al. (1999) when v = 0, i.e., for 
observations in a homogeneous universe. Then we employ the formulae presented here to obtain 
confidence contours in the £l m —v plane, constraining Q,^ so that = 1.0 — Vt m . By proceeding 
this way we are assuming that inhomogeneous matter, e.g., galaxies, are sufficiently removed from 
the lines of site of the 60 SNe and that lensing is negligible. 

Figure 1 presents the 68, 90 and 99% confidence contours for the fit (solid, long-dashed and 
short-dashed lines respectively). The results at v = are clearly consistent with the Perlmutter et 
al. (1999) findings. Unfortunately, no value of v can be ruled out with this sample because of its 
size and depth in redshift space. The best fit overall is at (Q m = 0.31, v = 0.0). 

A simple-minded way to estimate the number of additional SNe required to rule out any value 
of v at 99% confidence is to amplify the contribution of each SN to \ 2 by some factor. This has the 
effect of simulating more SNe which are exactly like the real ones, and hence the actual best fit will 
not move, but the contours will shrink. Figures 2 and 3 are the same contours as in Figure 1, but 
with sample size increased by factors of 10 and 100 respectively. We see that by simply enlarging 
the CT+SCP sample 10 times will allow a v = 2 value to be ruled out at 99% confidence. A 
factor of 100 will allow a much larger range of v to be excluded. These results are quite promising 
because a sample size of SNe la in the thousands extending to even higher redshifts should be 
possible with a ground based SNe pencil beam survey (Wang (2000b)) or a satellite mission such 
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as the SuperNova Acceleration Probe (SNAP - http://snap.lbl.gov). 

The fact that the current data at z < 1 do not rule out any value of v is not surprising. 
Assuming (fi m = 0.3, Oa = 0.7), the increase in distance modulus incurred by increasing v from 
to 2 at z = 0.5 is only 0.04 magnitudes (see Kantowski (1998)). However, at redshifts of z = 1.0, 
1.5 and 1.7 to be achieved in the future by SNAP, the increases are 0.14, 0.27 and 0.32 magnitudes 
respectively. Lensing complications will begin to occur at these higher redshifts and will have to be 
corrected for and/or selected against {see Wambsganss et al. (1997), Premadi et al. (1998), Holz 
& Wald (1998), Metcalf & Silk (1999), Tomita et al. (1999), and Wang (1999)}. Such data are 
likely to provide limited leverage in determining (fi m , Oa), unless v is properly constrained. Even 
though the most prevalent opinion is that v = 2, i.e., there is no significant intergalactic medium, 
there may be a sea of massive neutrinos etc., out there that makes < v < 2. 

As observations reach z ~ 3 lensing by galaxies becomes even more important and details 
of mass inhomogeneities will significantly distort the Hubble curve. Distances given here are still 
useful even without exact knowledge of those inhomoheneities. A lower bound on the primary im- 
age magnification at a given redshift (relative to the mean) is given by n cu toff = Dj(tt m ,£lA, v = 
0; z)/Dg(Q m , Qa, v\ z). This bound represents sources which, by chance, are not lensed. This num- 
ber can be compared to histograms given in Wambsganss et al. (1998), Premadi et al. (2001), Holz 
& Wald (1998), and Bergstrom et al. (2000). 

4. Conclusions 

We have given useful forms for the luminosity distance in the currently relevant inhomogeneous 
Qq = 1 FLRW cosmologies. These cosmologies are all dynamically FLRW in the large but differ in 
how gravitating matter affects optical observations. A beam filling parameter v, < v < 2 allows 
the matter to vary from completely transparent and homogeneous to completely inhomogeneous 
and exterior to any observing beam. In order to determine the values of the cosmic parameters (Sl m , 
Qa) from Hubble curves at high redshift, the value of v must also be constrained. For fixed values 
of (Q m ,^A), increasing v from (totally homogeneous universe) to 2 (totally clumped) increases 
the distance moduli of points on the Hubble curve, especially at higher redshifts as pointed out in 
§3. When observational error is taken into account, the problem of using standard candles at high 
redshift while ignoring v to obtain S7a will become particularly confounding. 

The current sample of SNe la at z < 1 fail to constrain the value of the beam-filling parameter 
v. Samples 10 to 100 times larger than the current sample, and in the same redshift range, will 
constrain v. In order to unambiguously determined (Q, m ,tlA) from even higher redshift observations 
like those planned in the future, the distance-enhancing effect of v must be accounted for in the 
luminosity distance formulae. 

The authors wish to thank Y. Wang, G. Kalbfleisch, David Branch, and E. Baron for exceed- 
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Fig. 1. — The Cl m -i/ plane with 68, 90 and 99% confidence contours (solid, long-dashed and short- 
dashed, respectively) resulting from an attempt to constrain those parameters using the 60 SNe la 
from the Calan/Tololo Supernova Search + Supernova Cosmology Project sample Perlmutter et 
al. (1999). The data were first fit assuming v = (completely homogeneous universe) to recover 
a result consistent with the original findings. Then assuming U m = 1 — J7a the above x 2 grid was 
calculated and contours of equal A% 2 above the minimum \ 2 were plotted. 

Fig. 2. — Same as Fig. 1 except 10 times more SN measurements exactly like the 60 SNe la from 
the Calan/Tololo + Supernova Cosmology Project sample were used in computing the confidence 
levels. 

Fig. 3. — Same as Fig. 2 except 100 rather than 10 times more SN were assumed to exist. 
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